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Shape complexity and fractality of fracture surfaces of swelled isotactic polypropylene 

with supercritical carbon dioxide 

Wei-Xing Zhou/'Q Bin Li, 1 Tao Liu, 1 Gui-Ping Cao, 1 Ling Zhao, 1 and Wei-Kang Yuan 1 

1 State Key Laboratory of Chemical Reaction Engineering, 
East China University of Science and Technology, Shanghai 200237, China 
(Dated: February 2, 2008) 

We have investigated the fractal characteristics and shape complexity of the fracture surfaces of 
swelled isotactic polypropylene Y1600 in supercritical carbon dioxide fluid through the consideration 
of the statistics of the islands in binary SEM images. The distributions of area A, perimeter L, 
and shape complexity C follow power laws p(A) ~ A _ ' ,1A+1 ', p(L) ~ L~' ML+1 - ) , and p(C) ~ 
C- (v+1 \ with the scaling ranges spanning over two decades. The perimeter and shape complexity 
scale respectively as L ~ A D ^ 2 and C ~ A q in two scaling regions delimited by A w 10 s . The 
fractal dimension and shape complexity increase when the temperature decreases. In addition, the 
relationships among different power- law scaling exponents fiA, /is, v, D, and q have been derived 
analytically, assuming that A, L, and C follow power-law distributions. 

PACS numbers: 05.45.Df, 61.25.Hq, 68.47.Mn, 82.75.-z 



INTRODUCTION 



Polypropylene, being one of the fastest growing en- 
gineering plastics, has wide industrial and everyday life 
applications due to its intrinsic properties such as low 
density, high melting point, high tensile strength, rigid- 
ity, stress crack resistance, abrasion resistance, low creep 
and a surface which is highly resistance to chemical at- 
tack, and so on 1]. Impregnation of nucleating agents 
allows the polymer to be crystallized at a higher temper- 
ature during processing and changes a lot the mechanical 
and optical performance of the resultant polymer matrix 
@i01 • In such nucleating agent impregnating processes, 
supercritical fluids are well-established swelling solvents 
whose strength can be tuned continuously from gas-like 
to liquid-like by manipulating the temperature and pres- 
sure. Especially, supercritical carbon dioxide is a good 
swelling agent for mostpolymers and will dissolve many 
small molecules QJa, |y, Q] , except for some fluoropoly- 
mers and silicones [8j ■ This provides the ability to control 
the degree of swelling in a polymer |jj Efl H3 as weu as 
the partitioning of small molecule penetrants between a 
swollen polymer phase and the fluid phase [LIL Il3| . 



The morphology of fracture surfaces of a material is 
of great concerns and interests in many studies |l4l 1 1 ol| . 
A frequently used tools is fract al g eometry pioneered by 
Mandelbrot's celebrated work [ly. Specifically, fractal 
geometry has been widely applied to the to pog raphical 
description of fracture surfaces of metals [r7tll8l . ceram- 
ics [H H El polymers |H H| ^oncretes 0, El IS 
E3, alloys Ha EllM H3 , rocks jH |H |H , and many 
others. There are many different approaches adopted in 
the determination of fractal dimensions of fracture sur- 
faces and the estimated fractal dimensions from different 
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materials vary from case to case and are not universal 

H 

Concerning the fractal characterization of the fracture 
surfaces of polypropylene, it was found that the fractal 
dimension of the fracture surfaces is less than 2.12 [23j| . 
In this paper, we investigate the shape complexity and 
fractality of fracture surfaces of nucleating agents im- 
pregnated isotactic polypropylene (Y1600) swelled with 
supercritical carbon dioxide based on the area-perimeter 
relationship [l6l fill l3C| . We shall see that the fractal 
dimension of the fracture surface of supercritical dioxide 
swelled polypropylene is much higher than that without 
swelling and decreases with increasing temperature. 



II. EXPERIMENTAL 

The isotactic polypropylene (Y1600) powder we have 
used has an average diameter of 3 ~ 4 mm. Fifteen grams 
of isotactic polypropylene powder was melt in an oven at 
200 °C and then made into film with a thickness of 0.3 
mm by using a press machine with a pressure of 0.5 GN. 
The size of the standard film samples was 1mm x 3mm 
and are refluxed with acetone for twenty four hours to 
remove the impurities, and then annealed at 200 °C for 
two hours. The nucleating agent we used was NA21. 

The schematic flow chart of the experimental appara- 
tus is illustrated in Fig. ^ The experimental apparatus 
consists mainly of a gas cylinder, a gas booster, a digital 
pressure gauge, an electrical heating bath, and valves and 
fittings of different kinds. The system was cleaned thor- 
oughly using suitable solvents and dried under vacuum. 
Isotactic PP films marked from one to twelve were placed 
in the high pressure cell together with desired amount 
(1% wt) of nucleating agent (NA21). The system was 
purged with CO2 and after the system had reached the 
desired temperature and thermal equilibrium, CO2 was 
charged until the desired pressure was reached. The im- 
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FIG. 1: Schematic flow chart of the experimental apparatus. 



pregnation process lasted for four hours and then de- 
pressurized the CO2 from the high pressure cell rapidly. 
Then the vessel was cooled and opened, and the spec- 
imens were taken out for analysis. All of the metallic 
parts in contact with the studied chemicals were made of 
stainless steel. The apparatus was tested up to 35 MPa. 
The total volume of the system is 500 mL. 

The experiments were performed at fixed pressure of 
7.584 MPa and at four different temperatures of T = 
60 °C, 80 °C, 100 °C, and 120 °C. For each sample, we 
took ten SEM pictures at magnification of 5000 or 10000. 
The pictures have 256 grey levels. Due to the nature of 
area-perimeter approach, the results are irrelevant to the 
magnification |16| . A grey-level picture was transformed 
into black-and-white image for a given level set e accord- 
ing to the criterion that a pixel is black if its grey level 
is larger than 256e and is white otherwise. The resultant 
binary image has many black islands. A typical binary 
image of islands is shown in Fig. El I n our calculations, 
we have used eleven level sets from 0.45 to 0.95 spaced 
by 0.05 for each SEM picture. 



III. FRACTALITY AND SHAPE COMPLEXITY 
A. Rank-ordering statistics of area and perimeter 

In order to estimate the probability distribution of 
a physical variable empirically, several approaches are 
available. For a possible power-law distribution with fat 
tails, cumulative distribution or log-binning technique 
are usually adopted. A similar concept to the comple- 
mentary distribution, called rank-ordering statistics |37j , 
has the advantage of easy implementation, without loss 
of information, and less noisy. 

Consider n observations of variable x sampled from a 
distribution whose probability density is p(x). Then the 
complementary distribution is P(y > x) — p(y)dy. 



FIG. 2: Image of the recognized islands from a typical SEM 
picture of swelled isotactic polypropylene impregnated with 
NA21 at level set of 0.7. 

We sort the n observations in non-increasing order such 
that x\ > X2 > ■ ■ ■ > xr > ■ ■ ■ > x n , where R is the 
rank of the observation. It follows that nP(x > xr) is 
the expected number of observations larger than or equal 
to xr, that is, 

nP(x >x R ) = R. (1) 

If the probability density of variable x follows a power 
law that p(x) ~ x _< - 1+AI - 1 , then the complementary dis- 
tribution P{x) ~ x^^ . An intuitive relation between xr 
and R follows 

x R ~ R- 1 ^ . (2) 

A rigorous expression of © by calculating the most prob- 
able value of xr from the probability that the i?-th value 
equals to xr gives [37) 

when fiR 3> 1 or equivalently 1 <C R < N, we retrieve 
• A plot of In xr as a function of In R gives a straight 
line with slope — l//i with deviations for the first a few 
ranks if x is distributed according to a power law of ex- 
ponent /!. 

We note that the rank-ordering statistics is nothing 
but a simple generalization of Zipf 's law [la and 
has wide applications, such as in linguistics [39j, the dis- 
tribution of large earthquakes |4fl |. time-occurrences of 
extreme floods 41], to list a few. More generally, rank- 
ordering statistics can be applied to probability distri- 
butions other than power laws, such as exponential or 
stretched exponential distributions |42|. normal or log- 
normal distributions |37l ]. and so on. 

Figure illustrates the log-log plots of the rank- 
ordered areas A at different temperatures. There is 
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FIG. 3: Log- log plot of the rank-ordered areas A at different 
temperatures shown in the legend. The plots are translated 
vertically for clarity. The solid lines are linear fits to the 
data at 59 < R < 3890 for T = 120 °C, 10 < R < 1995 for 
T = 100 °G, 22 < R < 1000 for T — 80 °C, and 8 < R < 800 
for T = 60 °C, respectively. 



FIG. 4: Log-log plot of the rank-ordered perimeters L at 
different temperatures marked in the legend. The plots are 
translated vertically for clarity. The solid lines are linear fits 
to the data at 5 < R < 3890 for T = 120 °C, 24 < R < 3000 
for T = 100 °C, 16 < R < 2000 for T = 80 °C, and 
10 < R < 1500 for T = 60 °C, respectively. 



clear power-law dependence between A and its rank R 
with the scaling ranges spanning over about two decades. 
Least-squared linear fitting gives 1/ha = 0.80 ± 0.05, 
1/fiA = 0.63 ± 0.02, l/(i A = 0.61 ± 0.03, and l/fi A = 
0.85 ± 0.04 with decreasing temperatures, where the er- 
rors are estimated by the r.m.s. of the fit residuals. Us- 
ing a^ A = Pi/ha/ we obtain the power-law exponents 
fj, A = 1.24 ± 0.08, (j,a = 1-59 ± 0.06, fi A = 1.64 ± 0.09, 
and [i a = 1.18 ± 0.05 with decreasing temperatures. 

Figure 0] shows the log-log plots of the rank-ordered 
perimeters L at different temperatures. There are also 
clear power-law relations between L and its rank R whose 
scaling ranges spanning over about three decades. Least- 
squared linear fitting gives 1/Vl = 0.56 ± 0.02, = 
0.51 ± 0.01, l/m = 0.55 ± 0.03, and l/[i L = 0.71 ± 
0.03 with decreasing temperatures. We thus obtain the 
power-law exponents hl = 1.78 ±0.07, /jl = 1.95 ±0.05, 
Hl = 1-83 ± 0.12, and //£ = 1-41 ± 0.06 with decreasing 
temperatures. It is noteworthy that the scaling ranges 
of the rank-ordering of both A and P are well above 
two decades broader than most of other experiments with 
scaling ranges centered around 1.3 orders of mag nitude 
and spanning mainly between 0.5 and 2.0 |43ll44| . 

B. Scaling between area and perimeter 

If an island is fractal, its area A and perimeter L follow 
a simple relation, 

L X I D ~ A 1 / 2 , (4) 

where D is the fractal dimension describing the wiggliness 
of the of the perimeter. For D = 1, the perimeter of the 
island is smooth. For D = 2, the perimeter becomes more 



and more contorted to fill the plane. This area-perimeter 
relation J3J can also be applied to many islands when 
they are self-similar in the sense that the ratio of L X I D 
over A 1 / 2 is constant for different islands 0|. It has 
been used to investigate the geometry of satellite- and 
radar-determined cloud and rain regions by drawing A 
against L for different regions in a log- log plot [3(|. The 
slope of the straight line obtained by performing a linear 
regression of the data gives 2/D. 

Following the same procedure, we plotted for each tem- 
perature L against A for all islands identified at different 
level sets. The resultant scatter plot for a given tem- 
perature shows nice collapse to a straight line. The four 
slopes are D = 1.34±0.08 for T = 60 °C, D = 1.32±0.07 
for T = 80 °C, D — 1.31 ± 0.08 for T = 100 °C, and 
D = 1.28 ± 0.09 for T = 120 °C, respectively. The errors 
were estimated by the r.m.s. of the fit residuals. 

However, a closer investigation of the scatter plots 
shows that there are two scaling regions with a kink at 
around A = 1000. To have a better view angle, we adopt 
an averaging technique for both A and L. We insert 
?i — 2 points in the interval [min(A), max(^4)] resulting 
in min(A) = ao < <xi < &2 < • • • < a„ = max(A) so 
that Oi's are logarithmically spaced. We can identify all 
islands whose areas fall in the interval [ai,et; + i). The 
geometric means of L and A are calculated for these is- 
lands, denoted as (L) and (A). The calculated means 
(L) and (A) are plotted in Fig. [3J The two scaling re- 
gions are fitted respectively with two straight lines and 
we have D x = 1.42 ±0.06 for T = 60 °C, D x = 1.33±0.06 
for T = 80 °C, D x = 1.33 ± 0.07 for T = 100 °C, and 
Di = 1.30 ± 0.05 for T = 120 °C, in the first region 
that (A) < 1000 and D 2 = 1.96 ± 0.07 for T = 60 °C, 
D 2 = 1.89 ± 0.24 for T = 80 °C, D 2 = 2.15 ± 0.11 for 
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FIG. 5: Plots of (L) with respect to (^4) at different tempera- 
tures. The plots are translated vertically for clarity. We note 
that (L) — 4 when (A) — 1 in all the four plots. 



T = 100 °C, and D 2 = 2.07 ± 0.14 for T = 120 °C in the 
second region that (A) > 1000. 

It is very interesting to notice that L oc A for large ar- 
eas. Under the consumption that the islands are fractals, 
this relation can be interpreted that the perimeters of the 
large islands are so wiggly that they can fill the plane. 
However, there are simple models that can deny this con- 
sumption of self-similarity. Consider that we have a set 
of strip-like islands with length I and width w. We have 
L = 2(1 + w) and A — Iw. When w is fixed and I » w, 
it follows that L oc A. This simple model can indeed ex- 
plain the current situation. For cracked polypropylene, 
there are big strip-like ridges in the intersection whose 
surface is relatively smooth. For small islands, the frac- 
tal nature is more sound. 

According to the additive rule of codimensions of two 
intersecting independent sets, the fractal dimension of 
the surface is D s = D x + 1 [l]|. Therefore, D s = 2.42 ± 
0.06 for T = 60 °C, D s = 2.33 ± 0.06 for T = 80 °C, 
As = 2.33 ± 0.07 for T = 100 °C, and D s = 2.30 ± 0.05 
for T = 120 °C. We see that the fractal dimension D s 
increases with decreasing temperature. In other words, 
with the increase of temperature, the surface becomes 
smoother. 

Assuming that p(A) ~ A^^ A+1 ^ and p(L) — 
we can estimate the fractal dimension D in 
(@J according to the relation p(A)dA = p(L)dL such that 



D = 2/i a /ij l 



(5) 



The four estimated values of D using this relation © are 
1.67, 1.79, 1.63, and 1.40 with increasing temperature. 
The discrepancy from D\ is remarkable (15%, 25%, 18%, 
and 7%, accordingly). The source of this discrepancy is 
threefold. Firstly, the power-law distributions of A and 
L have cutoffs at both ends of small and large values 



so that the derivation of JSJ is not rigorous. Secondly, 
the determination of the scaling ranges may cause errors. 
Thirdly, the scaling ranges of the three power laws are 
not consistent with each other. 



C. Shape complexity 

To quantify the shape complexity of an irregular frac- 
tal object, besides the fractal dimension, there are other 
relevant measures related to the fractal nature [ZBl l46l| . 
For a (i-dimensional hypersphere, its surface and volume 
are related by 



S, 



d,sph 



k,i 



(d-l)/d 
' d,sph 



where 



kd = 



r 1 / 2 (i + rf/2) 



(6) 



(7) 



Since a hypersphere has the maximal enclosed volume 
among the objects with a given surface Sd, we have 



1 



(d-l)/d 



Then the following dimensionless ratio 



kdSd 



V. 



(d-l)/d 



(8) 



(9) 



d 



describes the irregularity or complexity of the investi- 
gated object 0- It follows that 1 < Cd < oo. The lim- 
its are reached at Cd = 1 for hyperspheres and Cd = oo 
for fractals. For real fractal objects, sometimes known 
as prefractals the scaling range is not infinite. The 
shape complexity Cd is thus finite. In the case of d = 2, 
we have 



KC = 



27t1/2^1/2 



< OO 



(10) 



This dimensionless measure of shape complexity has been 
applied to scalar field of concentration in turbulent jet at 
high Reynolds numbers 0, H|| . 

In Fig. El is shown the log-log plot of the rank-ordered 
shape complexity C at different temperatures. The 
power-law dependence C ~ R^ 1 ^ implies a power-law 
probability density of shape complexity, that is, 



p{C) - c- (l/+1) 



(11) 



We find that v = 3.97±0.11 for T = 60 °C, v = 4.61±0.17 
for T = 80 °C, v = 4.93 ± 0.20 for T = 100 °C, and 
v = 5.22 ± 0.17 for T = 120 °C. The mean logarithmic 
complexity (InC) = l/{v — 1) is calculated to be 
0.34, 0.28, 0.25, and 0.24 with increasing temperature. 
In other words, the shape complexity of the islands de- 
creases with increasing temperature. This is consistent 
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FIG. 6: Log-log plot of the rank-ordered shape complexity C 
at different temperatures. The plots are translated vertically 
for clarity. The solid lines are linear fits to the data at 40 < 
R < 3890 for T = 120 °C, 24 < R < 5000 for T = 100 °C, 
40 < R < 3000 for T = 80 °C, and 25 < R < 1500 for 
T = 60 °C, respectively. 



with the fact the D decreases with increasing tempera- 
ture. 

Following the procedure in Sec. lIIIBl we calculate the 
means, (A) and (C). The results are shown in Fig. [7| 
Again, we see two scaling regions guided by straight lines 
fitted to the data 



C ~A q 



(12) 



The slopes are respectively q x 
0.17 for T = 80 °C, q x = 0.17 
0.15 for T = 120 °C, in the fir 
and q 2 = 0.54 for T = 60 °C, 
q 2 = 0.44 for T = 100 °C, and 
in the second region that (A) 
that qi — (Di — l)/2 with i = 
cases, which is nothing but a 
combination of Q and l|10|) . 

Similarly we can derive the 
fiA, and v as follows 



= 0.21 for T = 60 °C, q x = 
for T = 100 °C, and q x = 
st region that (A) < 1000 
q 2 = 0.58 for T = 80 °C, 
q 2 = 0.48 for T = 120 °C 
> 1000. We have verified 
= 1,2 holds exactly for all 
direct consequence of the 



relationship among q, D, 



2n A /(D - 1) = n A /q 



(13) 



Again, the discrepancy between the fitted v values and 
those estimated indirectly from i|13fl is remarkable (38%, 
48%, 53%, 30%), which can also be resorted to similar 
reasoning. 



IV. CONCLUSION 



FIG. 7: Plots of (C) with respect to (^4) at different tempera- 
tures. The plots are translated vertically for clarity. We note 
that (C) = 2/tt 1/2 when (A) = 1 in all the four plots. 



and perimeter of the islands are found to follow power 
laws spanning over two decades of magnitude, via rank- 
ordering statistics. The well-established power-law scal- 
ing between area and perimeter shows the overall self- 
similarity among islands when A < 10 3 . This shows that 
the fracture surface is self-similar. For large islands, the 
fractal dimension is estimated to be close to 2 which can 
be explained by the fact that most of the large islands 
are strip-shaped. The fracture surface is rougher at low 
temperature with larger fractal dimension. 

We have also investigated the shape complexity of the 
fracture surfaces using a dimensionless measure C. The 
distribution of shape complexity is also found to fol- 
low a power law spanning over two decades of magni- 
tude. The shape complexity increases when the island is 
larger. There are two power-law scaling ranges between 
C and A delimited around A — 10 3 , corresponding to the 
two-regime area-perimeter relationship. The exponent q x 
serves as an inverse measure of overall shape complexity, 
which is observed to increases with temperature. This is 
consistent with the change of fractal dimension at differ- 
ent temperatures. 

Furthermore, the relationships among different power- 
law scaling exponents of the probability distributions 
(pA, MBi ano - v ), of area-perimeter relation (-D), and 
of complexity-area relation (q) have been derived ana- 
lytically. However, these relations hold only when the 
probability distributions of A, L, and C follow exactly 
power laws. In the present case, we observed remarkable 
discrepancy between numerical and analytical results. 



We have investigated the fractal characteristics of the 
fracture surfaces of swelled isotactic polypropylene Y1600 
impregnated with nucleating agent NA21 through con- 
sideration of the statistics of the islands in binary im- 
ages. At a given temperature, the distributions of area 
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